1 Washington Department of Fish & Wildlife

Correspondence: Bethany Craig <>

1 Introduction

Let’s begin with a paragraph that’s a more general introduction in terms of why monitoring steelhead in the Dungeness is important? Part of a listed DPS, etc.

Steelhead spawning ground surveys in the Dungeness River basin are inherently challenging due to springtime snow melt and rain events which can lead to high, turbid water and unsafe survey conditions. In most years it is not possible to survey for steelhead through the entirety of the spawning season, and in some years poor survey conditions prevent an adequate number of surveys from being completed to produce an estimate of escapement based on redd counts. SONAR may provide an alternative method for steelhead enumeration and run timing in a dynamic, turbid snow-melt system like the Dungeness watershed. Since 2018 the Washington Department of Fish and Wildlife (WDFW) has operated a stationary multi-beam SONAR unit in the lower Dungeness River to enumerate and gather run-timing information on winter steelhead (Onchorhynchus mykiss).

2 Methods

2.1 SONAR operation

In 2018 the SONAR unit was deployed at approximately at river mile (RM) 0.2, below the majority of steelhead spawning habitat, and ensonified an approximately 20 meter (m) wide run in the river (Figure 2.1). The SONAR unit was mounted to a pole mount and attached to a reinforced ladder, secured to the river bottom by rebar (Figure 2.2).

Fish frequently milled or held in front of this SONAR site, which made counting fish passage difficult, and in 2019, the SONAR site was moved upstream to approximately RM 0.3, to a site with higher velocity, past which fish actively migrated (Figure 2.3). This site was easily accessible from the field trailer site, which enabled the unit to be directly connected and powered by trailer power, and any adjustments to the SONAR settings to be accomplished in the dry, safe comfort of the trailer.

Location of the SONAR site in the lower Dungeness River in 2018 (white) and 2019-2022 (striped).

Figure 2.1: Location of the SONAR site in the lower Dungeness River in 2018 (white) and 2019-2022 (striped).

SONAR unit deployment via a pole mount and ladder system in the Dungeness River in 2018.

Figure 2.2: SONAR unit deployment via a pole mount and ladder system in the Dungeness River in 2018.

Midway through the 2019 season the SONAR unit was mounted on a semi-permanent platform along the hardened left bank, in a spot that is protected and retains adequate depth so that the SONAR unit did not need to be shifted laterally to accommodate changing water levels (Figure 2.3).

In all years a picket weir was constructed approximately 1 m above upstream of the SONAR unit from the bank, extending out to approximately 1 m past the SONAR, to deflect debris (Figure 2.3). A second picket weir was constructed approximately 1 m below downstream of the SONAR unit to direct migrating fish out in front of the unit.

SONAR unit deployment via a pole mount and platform in the Dungeness River in 2020. Picket weir is the upstream picket weir.

Figure 2.3: SONAR unit deployment via a pole mount and platform in the Dungeness River in 2020. Picket weir is the upstream picket weir.

We deployed the ARIS 1800 Explorer, manufactured by Sound Metrics, of Bellevue, Washington. The ARIS 1800 uses 96 beams at 1.1/1.8 megahertz (Mhz) to project a 28-degree acoustic wedge. The SONAR unit was adjusted to have a pitch of 3.5 degrees to -8 degrees to ensonify the entire water column and channel. The unit was checked daily and adjusted as necessary to maintain full ensonification of the channel. Imagery was continuously recorded 24 hours a day, and saved in 30-minute files, so that 48 individual files were recorded for each full day of operation.

The SONAR was operated from early February or early March through late June or mid- to late July (Table 2.1, Figure 2.4). Other than a 22-day suspension in 2020 due to COVID-19 protocols, there were few outages and gaps in data collection (Table 2.1, Figure 2.4). SONAR imagery was reviewed for steelhead passage from the first day of operation through May 31th in each year.

Table 2.1: SONAR operational dates on the Dungeness River, 2018 - 2022.
Year First Date Last Date
2018 Mar 07 Jul 03
2019 Mar 05 Jul 17
2020 Feb 13 Jul 26
2021 Feb 01 Jun 22
2022 Feb 09 Jun 25
Plot of when SONAR was operational (purple) and not (yellow), from Feb 1 through July 31 each year.

Figure 2.4: Plot of when SONAR was operational (purple) and not (yellow), from Feb 1 through July 31 each year.

Plot of when SONAR was reviewed (purple) and not (yellow), from Feb 1 through July 31 each year.

Figure 2.5: Plot of when SONAR was reviewed (purple) and not (yellow), from Feb 1 through July 31 each year.

2.2 Data Processing

Because of the immense quantity of SONAR files, and the amount of time it took to process and review each 30-minute file, we initiated a subsampling scheme to enable the project team to complete review of the entire period of steelhead passage. In 2019 - 2022 the first 30 minutes of each hour were processed and reviewed for fish migration. A subset of days was fully reviewed (60 minutes of each hour) to compare fish migration with subsampled data (first 30 minutes of each hour). Another subset of days was double, or triple, reviewed to compare fish counts and lengths among each year’s 2 or 3 data reviewers. Table 2.2 shows what percentage of the hours were reviewed each year, and how much of the steelhead migration season was captured by operational SONAR.

Table 2.2: Number of hours during winter steelhead migration period (Feb 1 - May 31) and what percentage of those hours the SONAR was operational and what percentage of those hours have been reviewed.
Year Hours Operational (%) Reviewed (%)
2018 4,343 21.9 21.9
2019 4,343 71.1 40.0
2020 4,367 63.8 36.4
2021 4,343 73.4 41.9
2022 4,343 71.0 40.5

Each reviewed imagery file was processed using Sound Metric’s proprietary software ARISFish (v2.6.3 – v2.8.0). First, raw image files were background subtracted, which removed static objects from the image so that only objects in motion are shown. Then, an echogram was created, which transformed the image into a graph of distance (y-axis) and time (x-axis), so that objects in motion appeared as white “tracks.” The echogram enabled the data reviewer to quickly navigate to parts of the image file that contained objects that could be migrating fish. These tracks were then manually viewed alongside the raw image file to determine if the object was a fish to be further investigated.

Fish greater or equal to 45 centimeters (cm) were measured, marked, and counted using the ARISFish software. Forty-five cm was determined to be the minimum length of a potential steelhead, based on captures of steelhead during sampling in the Dungeness River 2014, 2015, and 2017 by the Jamestown S’Klallam Tribe (JSK) (unpublished data, C. Burns). Only fish that completely moved through the SONAR beams were counted; fish that nosed in and out or did not completely move from one side of the beams to the other were not counted. For each fish counted the following data were recorded:

  • Date
  • Hour of the 30-minute image file (e.g., 14:00, 14:30)
  • Time
  • Frame
  • Direction of travel (upstream or downstream)
  • Range (distance from the SONAR)
  • Length of the fish in cm
  • Data reviewer confidence (1 = extremely confident that the object counted is a fish \(\geq\) 45 cm, 2 = somewhat confident that the object is a fish \(\geq\) 45 cm, 3 = object of interest)

If no fish were observed in the 30-minute image file, a line of data with “NO FISH” was recorded to indicate that the file was reviewed for fish, but no fish \(\geq\) 45 cm were present. Marked fish were automatically saved within the image file for later error checking; data were also recorded within an Excel spreadsheet for data summarization and analysis.

2.3 Data Reviewer Comparison

In several years (2019 - 2021), a subset of SONAR footage was reviewed by all (two or three) of that years’ observers (Table 2.3). Within that subset of data, we summed the counts of each observer by date and direction, and examined the correlations between counts of different observers using the Pearson correlation coefficient. One observer (initial AS) was an observer every year, while other observers worked for one or two years only. We combined the daily counts across years when computing the correlation between observers. We also computed separate correlations for upstream and downstream counts.

To compare length measurements, we attempted to group individual observer detections of the same fish. We did this by first grouping fish detected in the same hour period moving in the same direction. If there was more than one fish detected in that group, we assumed that the relative lengths assigned by each observer corresponded to the same fish (i.e. the smallest observed fish by observer A was also the smallest observed fish by observer B). We then calculated the mean, standard deviation and coefficient of variation (CV) of the length measurements for each fish. Because observer AS was the most experienced, and measured nearly every fish, we also compared other observer length measurements against those by AS, and treated the measurements by AS as the benchmark for these comparisons. We summarized those differences with statistics like mean bias, mean absolute error (MAE), root mean squared error (RMSE) and mean absolute percent error (MAPE).

Table 2.3: Number of distinct days, hours and individual fish that were double or triple reviewed each year.
Year n Days n Hours n Fish
2019 5 32 54
2020 6 34 51
2021 12 65 92

2.4 Species Composition Sampling

Species composition sampling was conducted weekly, as river conditions allowed, during the period of SONAR operations in 2021 and 2022. A fine monofilament gill net 36 feet wide by 8 feet deep with a 2-inch mesh (4-inch stretch) was drifted through all sampleable habitat in the lower river from RM 3.3 to RM 0.5 (2021) or RM 0.8 (2022). Encountered fish were removed immediately from the net, sedated in a solution of tricaine mesylate (MS-222), and sampled for species, origin, length, gender, scales, and DNA. Captured steelhead were assessed for kelt status.

No regular species composition sampling was conducted in 2018, 2019, or 2020. In 2019, three sampling efforts targeting bull trout were conducted once per week in June at sites throughout the Dungeness and Gray Wolf rivers.

2.5 Determining Species

Bull trout are swimming past the sonar unit as well as steelhead, and we needed to parse which fish identified by the sonar were steelhead, and exclude any bull trout. The largest bull trout sampled in any species composition data was 67 cm, so we assumed any fish larger than 67 cm detected on the sonar was a steelhead. We then needed to determine of the fish equal to or less than 67 cm long, which of those are steelhead, and which are something else (e.g. bull trout, rainbow trout).

We only have species composition data for two years, 2021 and 2022. It was collected weekly, using tangle nets just upstream of the sonar location. For every fish caught, we know the date and fork length of that fish. Based on this data, we have also determined that the steelhead run on the Dungeness is over by early June. We made the assumption that any steelhead migrating up the Dungeness after June 1 would be balanced out by steelhead kelts that were observed moving downstream prior to that date. Therefore, we have only made predictions for fish detected prior to June 1st.

We modeled the probability of a fish being a steelhead using fork length and the Juilan day of capture. Since we only care about differentiating steelhead, we grouped resident rainbow and cutthroat trout with bull trout, and then fit a binomial GAM with a logit link, using splines of fork length and Julian day to predict the probability of a fish being a steelhead. We did not restrict the dataset to fish with fork lengths less than 67 cm, because larger fish have information about the shape of the logistic curve.

After fitting this GAM, we predicted the probability of being a steelhead for all fish observed on the sonar that were smaller than or equal to 67 cm, based on their length and Julian day of observation. Any fish with a probability of 50% or greater we assigned to be a steelhead.

2.6 Abundance Estimation

Although data was collected for the 2018 run year, the data quality was insufficient to estimate a steelhead abundance for that year. Frequent milling behavior at the SONAR site, counts that were much higher than subsequent years, and the fact that only 22% of the season was actually reviewed (Table 2.4) led us to treat 2018 as a pilot year and not make an escapement estimate.

For 2019 - 2022, SONAR fish targets that had a greater than 50% probability of being a steelhead, moved completely through the SONAR beams (direction of travel was upstream or downstream), and had a data reviewer confidence of 1 were included in the final steelhead counts and abundance estimate. We also restricted the steelhead season to SONAR observations or predictions between February 1 and May 31.

Table 2.4: Number of half hour periods each year between February and May, how many when the SONAR was operational, how many were reviewed and what percentage it was operational and reviewed, 2018 - 2022.
Year First Day Last Day n Periods n Operational n Reviewed Pct Operational Pct Reviewed
2018 Feb 01 May 31 5,758 1,470 1,465 25.5% 25.4%
2019 Feb 01 May 31 5,758 3,968 2,225 68.9% 38.6%
2020 Feb 01 May 31 5,806 3,566 2,331 61.4% 40.1%
2021 Feb 01 May 31 5,758 5,434 3,122 94.4% 54.2%
2022 Feb 01 May 31 5,758 5,020 2,870 87.2% 49.8%

To fill in the missing data from periods when the SONAR was not operating or the data had not been reviewed, this dataset was fit with a negative binomial generalized additive model (GAM) with a log link function using hour of day, discharge (from a USGS gauge on the Dungeness River, summarized at the hour time-scale) and Julian day-of-year as covariates. The GAM include a cubic spline for the hour of day, a thin-plate spline for discharge, and a factor spline for day of year, including an interaction between year and day of year. This GAM was then used to predict the upstream and downstream numbers of steelhead during all missing data periods, from February 1 through May 31 each year. This included all half hours when the SONAR was not operational for the entire 30 min, as well as all half hours that were not reviewed.

3 Results

3.1 Data Reviewer Comparison

After summing the counts across days when each pair of observers reviewed the same footage, Table 3.1 shows the total counts by year and direction (upstream vs. downstream) for each observer. Generally the upstream counts are very well aligned, and the downstream counts are close although they show more observer-to-observer variability than upstream counts.

Table 3.1: Total counts of steelhead by direction and year for each observer pair during days when both observers were counting.
Year Direction AS BC BT CS JG
2019 Upstream 38 38
2019 Upstream
38
37
2019 Upstream 38
37
2020 Upstream 31
32
2021 Upstream 14
13
2021 Upstream 26
26
2021 Upstream
3 3
2019 Downstream 6 7
2019 Downstream
7
5
2019 Downstream 9
7
2020 Downstream 15
12
2021 Downstream 19
19
2021 Downstream 24
20
2021 Downstream
6 1

The correlation coefficient between each pair of observers, by direction, is show in Table 3.2. Correlations for upstream moving fish were all greater than 0.97, while the correlation coefficients for downstream moving fish ranged from 0.17 to 0.87. Several of the weakest correlations did correspond with lower numbers of observed fish.

Table 3.2: Correlation (r) between daily counts of upstream and downstream moving fish for periods that were reviewed by pairs of observers. The total number of days when each pair reviewed SONAR data is shown, as well as the total number of fish observed by each observer.
Direction Obs 1 Obs 2 N Days Obs 1 Fish Obs 2 Fish r
Upstream AS BC 5 38 38 0.991
Upstream AS JG 5 38 37 0.998
Upstream AS BT 11 45 45 0.976
Upstream AS CS 6 26 26 1.000
Upstream BC JG 5 38 37 0.995
Downstream AS BC 3 6 7 0.866
Downstream AS JG 4 9 7 0.174
Downstream AS BT 10 34 31 0.863
Downstream AS CS 7 24 21 0.396
Downstream BC JG 3 7 7 0.500

For fish observed by multiple observers, the CVs of those multiple length measurements are displayed in Figure 3.1. The mean CV of length measurements for upstream moving fish was 0.06, while it was 0.09 for downstream moving fish. 95% of the upstream CVs were smaller than 0.15 (considered a reference point for precise measurements), while 84% of the downstream CVs met this criteria.

Coefficient of variation (CV) of length measurements for individual fish across multiple observers, colored by year and faceted by direction of movement. Dashed lines indicate the mean CV for that group, while the solid black line is a reference of 0.15.

Figure 3.1: Coefficient of variation (CV) of length measurements for individual fish across multiple observers, colored by year and faceted by direction of movement. Dashed lines indicate the mean CV for that group, while the solid black line is a reference of 0.15.

Histograms of length measurement bias, relative to observer AS, are shown in Figure ??. Summary statistics of the same comparison are depicted in Table @ref(tab:lgth_comp-tab). All observers had larger variability in their measurements of downstream moving fish compared to observer AS, as seen in RMSE, MAE and MAPE values. One observer (BT) had a positive mean bias relative to observer AS for both upstream and downstream fish, meaning they often measured the same fish as being slightly larger than AS. The other observers had negative mean bias, meaning they measured the same fish as being smaller compared to observer AS. This was more pronounced for downstream moving fish, up to an mean difference of nearly 12.5 cm per fish for observer BC (although that observer only measured four downstream fish).

The mean absolute percent error (MAPE) ranged from 6.7 - 9.3% for upstream moving fish and from 11 - 14.5% for downstream moving fish. To put that on the scale of fish lengths (centimeters), the root mean square error (RMSE) ranged from 6.7 - 8.3 cm for upstream moving fish and 8.5 - 17 cm for downstream moving fish.

Histograms and density plots showing the bias of length measurements, relative to observer AS, colored by observer and faceted by direction and year.

Figure 3.2: Histograms and density plots showing the bias of length measurements, relative to observer AS, colored by observer and faceted by direction and year.

Table 3.3: Summary statistics comparing length measurements of different observers with those made by observer AS.
Direction Observer N Fish Mean Bias RMSE MAE MAPE
Upstream BT 41 1.5 8.3 5.7 9.1
Upstream JG 37 -6.3 8.3 6.8 9.3
Upstream BC 36 -2.3 6.7 4.9 6.7
Upstream CS 25 -2.4 6.7 5.2 7.3
Downstream BT 24 1.5 9.7 7.4 12.6
Downstream CS 12 -7.5 8.5 8.0 11.0
Downstream JG 5 -7.4 12.1 10.7 13.0
Downstream BC 4 -12.4 17.0 13.0 14.5

3.2 Species Composition

Thirty-three species composition sampling events were conducted between early February and late June across 2021 and 2022. Steelhead and bull trout were encountered during sampling. The majority of captured steelhead were natural-origin; three hatchery-origin steelhead were captured in early April 2021, and four hatchery-origin steelhead were captured in 2022- three in mid- March and one in late May. Steelhead kelts were encountered starting in May 2022, and in early June 2021 (Table 3.4. From the species composition netting, there are 92 fish to use in this model. The lengths of these fish can be seen in Figure 3.3.

Table 3.4: Counts of species encountered during species composition sampling in the lower Dungeness River in 2021 and 2022.
Year Date Bull Trout Steelhead Resident Rainbow Cutthroat
2021 Feb 17 1 0 0 0
2021 Feb 22 0 2 0 0
2021 Feb 25 0 2 0 0
2021 Feb 26 0 1 0 0
2021 Mar 03 0 1 0 0
2021 Mar 04 1 0 0 0
2021 Mar 17 0 1 0 0
2021 Mar 18 0 1 0 0
2021 Mar 25 0 2 0 0
2021 Mar 26 0 1 1 0
2021 Apr 01 1 21 0 0
2021 Apr 06 2 0 0 0
2021 Apr 07 0 2 0 0
2021 Apr 08 0 61 0 0
2021 Apr 14 0 4 0 0
2021 Apr 15 3 0 0 0
2021 Apr 16 1 1 0 0
2021 Apr 28 0 1 0 0
2021 May 05 0 1 0 0
2021 Jun 08 4 1a 2 0
2021 Total 13 29 3 0
2022 Feb 02 3 0 0 1
2022 Feb 15 1 0 0 0
2022 Feb 22 0 2 0 0
2022 Mar 09 1 22 0 0
2022 Mar 16 0 32 0 0
2022 Mar 23 1 1 0 0
2022 Mar 30 0 2 0 0
2022 Apr 07 1 3 0 0
2022 Apr 19 2 7 0 0
2022 May 04 2 6b 0 0
2022 May 20 2 32b 0 0
2022 Jun 01 2 1b 0 0
2022 Jun 22 1 0 0 0
2022 Total 16 30 0 1
1 1 hatchery-origin steelhead was captured on April 1, 2021; and 2 hatchery-origin steelhead were captured on April 8, 2021.
2 1 hatchery-origin steelhead was captured on March 9, 2022; 2 hatchery-origin steelhead were captured on March 16, 2022; and 1 hatchery-origin steelhead was captured on May 20, 2022.
a 1 steelhead was identified as a kelt on June 8, 2021.
b 4 steelhead were identified as kelts on May 4, 2022; 2 steelhead were identified as kelts of May 20, 2022; and 1 was identified as a kelt on June 1, 2022.
Histogram of forklengths, colored by species. Dashed lines are the overlap range between steelhead and other species.

Figure 3.3: Histogram of forklengths, colored by species. Dashed lines are the overlap range between steelhead and other species.

3.3 Determining Species

The impacts of Julian day and fork length on the probability of a fish being a steelhead are displayed in Figure 3.4. Generally, larger fish have a greater chance of being a steelhead, but the exact break point of length (i.e. 50% probability of being a steelhead) shifts over the course of the season.

Fish greater than 67 cm were deemed steelhead, regardless of the Julian day of year they were observed. Fish less than or equal to 67 cm were deemed steelhead if the predicted probability of being a steelhead from our GAM model was greater than or equal to 50%. The number of fish within each category is shown in Table 3.5.

Julian day and fork length of each fish detected on SONAR, colored by the predicted probability of being a steelhead. Size of the point corresponds to the precision of the predicted probability (larger points indicate greater precision). Dark line indicates the 50% probability of being a steelhead.

Figure 3.4: Julian day and fork length of each fish detected on SONAR, colored by the predicted probability of being a steelhead. Size of the point corresponds to the precision of the predicted probability (larger points indicate greater precision). Dark line indicates the 50% probability of being a steelhead.

Table 3.5: Number of fish observed moving downstream or upstream (n Fish), less than or greater than 67 cm, and whether they were likely (> 50% probability) to be a steelhead. Totals within each movement and size group are shown (n Sthd) as well as totals for each year by direction (Total Sthd).
<= 67 cm
> 67 cm
Direction Year n Fish Likely Sthd n Sthd n Fish Likely Sthd n Sthd Total Sthd
Downstream 2019 79 49 49 180 180 180 229
Downstream 2020 226 144 144 119 119 119 263
Downstream 2021 186 108 108 213 209 213 321
Downstream 2022 188 119 119 105 104 105 224
Upstream 2019 206 131 131 459 459 459 590
Upstream 2020 411 279 279 406 406 406 685
Upstream 2021 243 148 148 492 469 492 640
Upstream 2022 438 292 292 269 267 269 561

3.4 Abundance Estimation

The estimated marginal effects of various covariates on the expected number of upstream and downstream moving fish are shown in Figures 3.5 and 3.6 respectively. For both directions, the model indicates that a greater number of steelhead will move during the night, especially in the late evening, compared to the daytime, and more will move when discharge flows are near 600 cfs. There is also a clear effect of day-of-year, with upstream numbers peaking near the beginning of April, while downstream numbers peak near the beginning of May, with some year-to-year variability in that run-timing curve.

Estimates of annual upstream and downstream migration of winter steelhead are shown in Table 3.6, as well as estimates of total net upstream escapement.

Table 3.6: Annual estimates of winter steelhead upstream, downstream, and net escapement to the Dungeness River between Feb 1 and May 31.
Upstream
Downstream
Net Escp.
Year Estimate SE Low CI Upp CI Estimate SE Low CI Upp CI Estimate SE Low CI Upp CI
2019 1,297 246.3 1,144 2,013 536 240.7 443 1,336 748 72.3 616 921
2020 1,633 105.9 1,460 1,868 669 47.3 593 777 963 68.4 846 1,112
2021 1,179 27.1 1,131 1,234 606 19.3 572 648 571 12.2 550 596
2022 1,059 30.6 1,009 1,123 436 22.6 407 488 622 13.1 599 648
Marginal effects for upstream-moving fish as estimated by the GAM for hour of day (A), discharge (B), and Julian day of year (C). Shaded areas indicated 95% confidence intervals.

Figure 3.5: Marginal effects for upstream-moving fish as estimated by the GAM for hour of day (A), discharge (B), and Julian day of year (C). Shaded areas indicated 95% confidence intervals.

Marginal effects for downstream-moving fish as estimated by the GAM for hour of day (A), discharge (B), and Julian day of year (C). Shaded areas indicated 95% confidence intervals.

Figure 3.6: Marginal effects for downstream-moving fish as estimated by the GAM for hour of day (A), discharge (B), and Julian day of year (C). Shaded areas indicated 95% confidence intervals.

4 Discussion

We used SONAR to make annual estimates of winter steelhead on the Dungeness River for spawn years 2019 - 2022. To do so, we needed to collect and utilize independent species composition data, and develop an analysis technique that allowed us to estimate passage during periods when the SONAR was not operable or the video had not been reviewed. In addition, analyzing these initial years of data has prompted several questions worth mentioning.

First, the variability between observers was relatively small. The total counts between pairs of observers each year were very similar, especially for upstream moving fish (Table 3.1). The coefficient of variation among different observers’ length measurements was quite small, suggesting good consistency across observers. It should be noted that we matched up which fish each observer measured in such a way as to minimize observer-to-observer differences for multiple fish observed in the same half-hour period. The overall bias in measurements, compared to observer AS, was quite low, but not zero. This could have an impact on whether a fish is predicted to be a steelhead or not, since length is a determining covariate in the species prediction model. However, Figure 3.2 suggests that there are differences in length measurements in both the positive and negative direction, so the impact on steelhead counts should be minimal.

Second, SONAR location is a crucial component of a successful SONAR project. After observing the fish behavior at the initial location in 2018, the decision was made to move the SONAR to a habitat unit that was more of a riffle, conducive to fish moving through rather than milling around.

Bethany, anything to add here about location choices?

In addition to the location, several factors in 2018 led to unrealistic estimates of escapement. The amount of time when the SONAR was either not in operation or that went un-reviewed in 2018 was much higher than in other years (Figures 2.4 and 2.5 and Table 2.4). In addition, the counts of both upstream and downstream moving fish were much higher in 2018, presumably due in part to more milling behavior at that site. The mean count of upstream moving fish in 2018 in a half-hour period was 1.02, while it ranged from 0.24 - 0.35 in 2019 - 2022. For downstream moving fish it was a mean of 0.70 in 2018, and a range of 0.1 - 0.15 in 2019 - 2022. Due to these reasons, we made the decision to treat 2018 as a pilot year and not make an escapement estimate using SONAR that year, and instead utilize the lessons learned from that year to improve operations in subsequent years.

Third, how to best deal with steelhead kelts is an open question. For fish that move upstream, downstream, and upstream again before spawning, we would like to subtract the downstream counts from the upstream counts, to avoid double-counting the same fish. However, steelhead that kelt may be observed moving upstream and again moving downstream, and we would not like those counts to cancel each other out. Currently, we are sidestepping this issue by assuming that the number of kelts detected prior to June 1 are balanced out by upstream moving steelhead after that date. However, that date is at best an educated guess which could lead to bias in our estimates in one direction or the other. Continued investigation of this issue is warranted, perhaps by more closely examining the species composition data for kelting behavior (or changing that sampling to better trap kelts) or by analyzing metrics such as the upstream:downstream count ratios through the season, to determine if there is a pattern that could guide what that equilibrium date should be.

Finally, the species composition data, while an extremely important part of this SONAR project, did raise a few issues for future consideration. A few bull trout were captured multiple times in the same pool during species composition netting, suggesting that the species composition sampling may be sampling not just fish that are moving, but also fish that may be holding (i.e. not moving past the SONAR). Whether this is actually impacting the predictions of species based on size and Julian day is unclear. In addition, the current methods cannot account for potential differences in what hour of the day certain species are more or less likely to move. If steelhead are more likely to move at night, we would want to incorporate hour of day into our species prediction. Without a different study design for species composition, we are unable to address this question. However, it seems reasonable to assume that steelhead and bull trout have similar behavior with regards to movement and time of day.

Should we bring anything into the discussion about comparisons between SONAR estimates and redd estimates of escapement?

4.1 Recommendations

We recommend in the future trying to set the SONAR unit up earlier in the year, before the steelhead start moving into the Dungeness. Even in 2021, when the SONAR was deployed starting February 1, there was still one steelhead detected on that date, suggesting the run may begin sooner than that. If the SONAR is deployed earlier, then we would also recommend that the species composition sampling begin at the same point in the year that the SONAR is deployed, to ensure that we can identify when steelhead (as opposed to other species) are moving in the Dungeness. After a few years of deploying quite early, we may begin to move that deployment date forward with a firmer understanding of when the run starts.

Another recommendation is to investigate the use of image recognition software and machine learning processing to help automate the review of SONAR images. Technicians are required to spend a substantial amount of time reviewing SONAR video, and if that could be reduced it would have multiple benefits. First, cutting hours of tedious review time could potentially be cost-effective, and make the data available for analysis sooner. Second, it would make the counts more reproducible and less subject to an individual observer’s skill and experience. Finally, by automating the review process there would no longer be a need to selectively choose which periods to review (e.g. first 30 minutes of every hour). Instead, counts could be obtained for the entire period the SONAR is deployed, reducing the reliance on estimates for periods with missing data.

5 Acknowledgements

Andrew Simmons, Dan Gorze, Jayson Gallatin, Brent Trim, and Greg Boecker of WDFW operated, maintained, and reviewed SONAR imagery data. Calvin Stokes also reviewed SONAR imagery. Peter Topping provided project support and advice. Kathryn Sutton led weekly species composition sampling with the help of many volunteers including Calvin Stokes, Jenni Whitney, Scott Williams, Jeff Gufler, David Thomas, Jacob Portney, Caleb Owen, and Kayla Own of WDFW, Chris Burns, Jarrett Burns, Casey Allen, and Aaron Brooks of the Jamestown S’Klallam Tribe, Dave Shreffler of Shreffler and Associates, Keith Denton of Keith Denton and Associates, and Roger Peters, Roger Tabor, Greg Byford, Evan Lewis, Jakob Bengelink, and Jack Brill of the U.S. Fish and Wildlife Service.